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Spatial domain enhancements 


= Transformations 


o intensity transformations 
° negatives 
° logs 
° power-law (gamma) 
e contrast stretching 
e level slicing 


e bit-plane slicing 


o distribution transformations 


e histogram equalization 


= Spatial filtering 
o image filtering 
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Logs 


u Log transformations 
o used to expand values of dark pixels q = 
o simultaneously compressing bright pixels a — pra kk, 
o compresses dynamic range of images A 


e Fourier spectrum 
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Gammas 


= Power-law transformations 


o sensors respond according to power law 
e CMOS, scanners, printing, displays 


e CRT: intensity to voltage response as power function (y’ = 1.8~2.5) 


o gamma correction 
o device dependent y 
o y variation also varies the color ratios 


o correct color reproduction needs knowledge of y 


o gamma injection 


o post image processing for contrast manipulation 
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Gammas 


= Power-law transformations 


L-1 —— EI 
o sensors respond according to power law ert 
e CMOS, scanners, printing, displays y= 0.10 
e CRT: intensity to voltage response as power function (y’ = 1.8~2.5) 3L/4 y = 020 
d y = 0.40 
o gamma correction > ME A 
o device dependent v E L/2 ġel 
o y variation also varies the color ratios š y = 1.5 
o correct color reproduction needs knowledge of y 6 y = 25 
LA y = 50 
o gamma injection y = 10.0 
o post image processing for contrast manipulation y = 25.0 
% L/4 IA 3L/4 L-1 


Input intensity level, r 
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Gammas eez? 


y correction 


j 


WITHOUT GAMMA CORRECTION WITH GAMMA CORRECTION 


credit: blog.wolfire.com 
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Gammas 


= y injection 
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Transformations 


= Compositions 
o piecewise combinations 


o piecewise linear 
e many T; formulated with this 


° need more user input paras 
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Output intensity level, s 


3L/4 


LO 


L/4 


Negative 


L/4 


Inverse log 


L/2 


Input intensitv level, r 


3L/4 


L-1 
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Contrast stretching 


= Contrast 


o low contrast images 
° due to poor illumination, low dynamic range sensors 


e wrong setting of lens aperture 


o full range stretching 
° (ms) = (Tmin: 0) 


° (15,52) = (maxi — 1) 


o thresholding 
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Contrast stretching 


= Contrast 


o low contrast images 
° due to poor illumination, low dynamic range sensors 


e wrong setting of lens aperture 


o full range stretching 
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pixel ranges 
° (%,S2) = (max L — 1) 


o thresholding 


EE604: IMAGE PROCESSING sandhanQiitk.ac.in 


Contrast stretching 
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o low contrast images 
° due to poor illumination, low dynamic range sensors 
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o full range stretching L-41 
kl Y 
i as 2 pixel ranges 
° (2,52) = (Tnax L — 1) 8 3L/4 
> 
£ 1/2 
E 
o thresholdin E 
B > L 
° = 
° 
0 


0 LA L/2 3L/4 L-1 


Input intensity level, r 
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Contrast stretching 


= Contrast 


o low contrast images 
° due to poor illumination, low dynamic range sensors 


e wrong setting of lens aperture 


o full range stretching L-41 
j (ri, 51) = (min 0) be f š 
> ixel ranges 
° (72,52) = (max: L — 1) Š 3L/4 p g 
E 
3 Eje SEM image of pollen grains 
E sa 
hresholdin $ 
o thresholding B La 
° Ti =n, Saz JU Saz L — 1 a 
0 


0 LA L/2 3L/4 L-1 


Input intensity level, r 
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Level slicing 


= Intensity levels 
o local thresholding, stretching 
o enhancing only specific intensities 


° e.g. detecting water, wetland in sat. images 
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Level slicing 


a Intensity levels 


o local thresholding, stretching 


o enhancing only specific intensities `L ST 


° e.g. detecting water, wetland in sat. images R 
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Bitplane slicing 


= Bitplanes 


o contribution of each bit for 
total image appearance 


o gives clue for a compression 
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One 8-bit byte pra 


Bit plane 8 
(most significant) 


Bit plane 1 
(least significant) 
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Bitplane slicing 


= B itp | anes SA a PA Bit plane 8 


(most significant) 
o contribution of each bit for 
total image appearance Bit plane 1 


(least significant) 


o gives clue for a compression 


o slicing 
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Bitplane slicing 
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(most significant) 
o contribution of each bit for 
total image appearance Bit plane 1 


(least significant) 


o gives clue for a compression 


o Slicing o reconstruction 


° bitplanes (8+7) 


bitplanes (8+7+6) 


bitplanes (8+7+6+5) 


EE604: IMAGE PROCESSING sandhanQiitk.ac.in 


Spatial domain enhancements 


= Transformations 


o intensity transformations g(x,y) = T (f(x, y)) 
* negatives bar egar 
° logs 
° power-law (gamma) ser 


e contrast stretching 
e level slicing 


e bit-plane slicing 
o distribution transformations g(x,y) =T; DS y))) 


e histogram equalization 


= Spatial filtering 
o image filtering 


EE604: IMAGE PROCESSING sandhanQiitk.ac.in 


Histograms 
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Histogram equalization 


= Assume fere La 1 d 
o T(r) is monotonic T s = T(r) — 
o bounded 0 < TG) <L —1 TO 
o variable equivalence sk 


° to cover all notations 


L-1 
T(r) 
4 
L=1|-=-=========== | 
Single | 
value, sy | 
TO i 
Single | 
value, s; | 
| 
| 
| 
| 
| 
| 
—— r 
0 Multiple Single L — 1 


values value 
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Y=T(X) TO) 
= Assume O<sr<L-1 
o T(r) is monotonic T S = T(r) det ps 
o bounded 0 < T(r) < L — 1 / \ Td — 
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Sk 
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u Assume QO=r=L-1 
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HAGO T(r) is cts & differentiable | 
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= cumulative function satisfies n 
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0 | 
0 Multiple single L-1 f 


values value 


EE604: IMAGE PROCESSING sandhanQiitk.ac.in 


Histogram equalization 


Y=T(X) T(r) 
u Assume OQerelL—i 
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HAGO T(r) is cts & differentiable | 
já É ps(s) pr) | 
e to cover all notations l 
py (y) px(x) | 
7 * 
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= cumulative function satisfies 4 
above properties for T(r) r | A 
s=T(r)=(L—1)] p(w) dw Single |. | 
0 value, sy | 
k TO 
sk = T(n) = (L — 1) > p,(rj) x Ge | 
j=0 Y= T(X) = (L = »/ px(z)dx i | 
0 | 
_ L - D< _ _ | 
= an R=EULA EI Melo : 
0 Multiple Single L — 1 
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Histogram equalization 


y LES T(X) is cts & differentiable 


What is py(y)? 


EE604: IMAGE PROCESSING 


Histogram equalization 


y LES T(X) is cts & differentiable 


What is py(y)? 


y 
/ py(z)dz = probability that 0 < Y < y 
0 


EE604: IMAGE PROCESSING 


Histogram equalization 


FETA) T(X) is cts & differentiable 


What is py(y)? 


y 
/ py(z)dz = probability that 0 < Y < y 
0 


= probability that 0 < X < T”*(y) 


EE604: IMAGE PROCESSING 


Histogram equalization 


x FETA) T(X) is cts & differentiable 
Y =T(X) = (L= »/ px (z)dz 
0 
What is py(y)? 


y 
/ py(z)dz = probability that 0 < Y < y 
0 


= probability that 0 < X < T”*(y) 


TTI 
= / px (w)dw 
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Histogram equalization 
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Histogram equalization 
Yaa] T(X) is cts & differentiable 


E =T(X)= L= DJ px (z)dz / N 


What is py(y)? 


i (f'm) = p(T") (TW) 


py (y) 
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Histogram equalization 


FETA) T(X) is cts & differentiable 


dl G i 
= Pn da PTW q E (y)) 
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Histogram equalization 


LEZA E T (X) is cts & differentiable 


dT do, 
= TET aero) 


wa 
ela 
a 
— 
b 
L 
— 
<= 
Na 
A 
| 
=| 
= 
<= 
| 
bek 
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Histogram equalization 


LEZA E T (X) is cts & differentiable 


dT to. 
méri da PTO GF '(9)) 2T(T(y)) = £y=1 
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Histogram equalization 


y LES T(X) is cts & differentiable 
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Histogram equalization 


Yaa] T(X) is cts & differentiable 


px(x) py (y) 
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Histogram equalization 


Histogram equalization 


Histogram equalization 
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Histogram equalization 


6 
< 
f 
Ea 
Lu 
O 
te) 
E 
(aL, 
Lu 
5 
< 
= 
E 
= 
© 
Lu 
Lu 


Histogram equalization 


Histogram equalization 


MW 77 $ 


LA DD L 


Histogram ezere, O r 
equalization 
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Histogram 


equalization 


= Color conversions 


o colors can be mapped 
with certain functions 


o Mapped images then 
histogram equalized 


EE604: IMAGE PROCESSING 


LEE töltő Color functions & Histogram equalization results 
E a EEE fa(RGB) 
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Histogram equalization 


Histogram equalization 
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Histogram equalization 


= Global 


Original iamge Histogram of original iamge 


10000 


0 100 200 300 


Image constructed using Equalized Histogram Equalized Histogram 
pr 
10000 


GE 


IK 


4 
Ore o, 
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liv: b 
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0 100 300 
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Histogram equalization 


= Global 


Original iamge Histogram of original iamge 
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= Local 
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Histogram equalization 


= Global 


Original iamge Histogram of original iamge 


10000 


0 100 200 300 


Image korsi Paċiea using Equalized Histogram Equalized Histogram 
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= Local 
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Histogram equalization 


Ref: wikipedia 
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Histogram equalization 


= Bilinear interpolation 


Ref: wikipedia 
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Ref: wikipedia 
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Histogram equalization 


= Bilinear interpolation 


Hen) S Nn) + E AQ), 


Ta — Z to — T 
f(x, y2) = — f(Qi2) + aa f(Q2). 


Ref: wikipedia 
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Histogram equalization 


= Bilinear interpolation 


f(z, m) = f(Qu) + — = f(Qa), 
flay) = 2 — HQu) gri H(Qu). 
BB (2 Qn) + EL (Qp) 


l «rrura ota — y) + eats egite 


` aaa erailak pl =" 


y) + f(Qiz)(r2 — 2)(y — vi) + f(Q22)(z giar y1)) 


EE604: IMAGE PROCESSING 


Histogram equalization: CLAHE 


Histogram equalization: CLAHE 


= ARE 
o Adaptive hist eq 


Histogram equalization: CLAHE 


= AHE =" CL 
o Adaptive hist eq o Clip limit 


Histogram equalization: CLAHE 


= AHE = CL = Interpolation 
o Adaptive hist eq o Clip limit o Bilinear 


Histogram equalization: CLAHE 


AHE CL Interpolation 
Adaptive hist eg Clip limit Bilinear 
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Histogram equalization: CLAHE 
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Conclusion 


Intensity transforms 


Distribution transforms 


Intensity transformations 


negatives 

logs 

power-law (gamma) 
contrast stretching 
level slicing 


bit-plane slicing 


Distribution transformations 


Histogram equalization 
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